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Abstract 

We investigate numerically a single-pulse implementation of a quantum Control- 
Not (CN) gate for an ensemble of Ising spin systems at room temperature. For an 
ensemble of four-spin "molecules" we simulate the time-evolution of the density 
matrix, for both digital and superpositional initial conditions. Our numerical 
calculations confirm the feasibility of implementation of quantum CN gate in 
this system at finite temperature, using electromagnetic vr-pulse. 
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Since 1994, Ising spin systems have been considered to be promising candidates for 
implementing quantum logic gates and for quantum computation It was shown 

in 0, ^ that quantum computation can be performed effectively using an ensemble of 
Ising spin systems at finite temperature. Recent implementations of two-qubit quantum 
logic gates in an ensemble of spin systems require application of complicated sequences 
of electomagnetic pulses to the nuclear spins with very close resonant frequencies 0, ||. 
We consider a simple single-pulse implementation of two-qubit quantum logic gate 
which can be realized in the systems with significant difference between the resonant 
frequencies. In order to understand the behavior of these systems, it is very important 
to perform numerical experiments of the dynamics of the single-pulse quantum logic 
gates and quantum computation for an ensemble of Ising spin systems at room tem- 
perature. Recent suggestions for quantum computation assume the operations using 
quantum superpositional (entangled) states. The process of manipulating these states 
in theoretical studies is different from the process of implementation of these states 
in real physical systems. There are two main problems. One problem is that these 
superpositional states are not eigenstates of the corresponding Hamiltonians. These 
superpositional states constantly change in time. One cannot create time-independent 
superpositional states even for an isolated system. The second problem is that the 
dynamics of quantum computation involves both resonant and non-resonant interac- 
tions. Non-resonant interactions are usually ignored, but they can inhibit the desired 
effects. For many physical problems (for reasonably small interactions), the resonant 
dynamics dominates. In this case, the non-resonant dynamics plays an insignificant 
role which can be reasonably estimated. The situation is completely different for quan- 
tum computations because non-resonant effects can accumulate with time, and create 
significant errors. 

In this paper, we present a numerical analysis of the quantum Control-Not (CN) 
gate at room temperature. We show that when the initial spin's frequencies differ signif- 
icantly and the constant of interaction between spins is small enough, the non-resonant 
effects give small contribution to the dynamics, allowing the single-pulse quantum CN 
gate to be implemented in this system. 

Quantum CN gates are of the central importance for quantum computation. Any 
quantum logic gate can be decomposed into a set of one-qubit rotations and CN gates 
[0]. The two-qubit gate can be described by the operator, 

CN = |00 >< 00| + |01 >< 01| + |10 >< 11| + 111 >< 10|. (1) 

The first number, i, in \ij > refers to the control qubit, which does not change its 
value during the CN operation. The second number, j, refers to the target qubit, 
which changes its value only if the control qubit has value "1". The simplest CN gate 
can be implemented using only one vr-pulse [§1-|10[- It was shown in |]T0[ that, in fact. 



a TT-pulse provides a modified CN gate, which differs from the "pure" CN gate (1) by 
a phase shift. Nevertheless, the desired feature of the CN gate remains: the vr-pulse 
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changes the state of the target qubit only if the control qubit has value "1". 

Consider an ensemble of four-spin molecules at finite temperature. Each molecule 
consists of four spins I = 1/2; each spin interacts with the other three spins, through 
an Ising interaction. The whole ensemble can be described by the density matrix, 
p. The Hamiltonian of this system, in the reference frame rotating with frequency uo, 
circularly polarized in xy plane magnetic field, h = h{cosujt, — sin cut), is. 



n 



3 
a=0 



In (2), 



0J„ 



b>a 



-fBa, ^ = 'yh. 



(2) 



(3) 



where 7 is the gyromagnetic ratio of each of the four spins; Ba is the permanent 
magnetic field which points in the positive 2;-direction, and is supposed to be different 
for different spins; Q is the Rabi frequency; Jab is the Ising interaction constant. We 
shall consider the case in which Jab = J- The energy levels for a single molecule 
containing four spins, for h = u = 0, are shown schematically in Fig. 1. The energy of 
the ground state is, 

^o = -^(E^a + 6j). (4) 

The first four excited states in Fig. 1 correspond to the one-spin excitations; the next 
six states correspond to the two-spin excitations; the next four states correspond to the 
three-spin excitations, and the last state, |1111 >, corresponds to the total inversion 
of the spin molecule. For the density matrix formalism, it is convenient to use the 
decimal notation. In this notation we have, 



|0000 >= IO3O2O1O0 >^ |0 >, |0001 >= IO3O2O1I0 >^ |1 >, 



115 >, 



(5) 



11111 >=|l3l2lllo>- , 

as shown in Fig. 1. Following the idea suggested in we assume that (using a 
particular sequence of electromagnetic pulses) an ensemble of four-spin molecules is 
initially prepared in the state which can be described by the density matrix. 



p = E/16 + pA, 



(6) 



Pa 



2kBT 



|0 >« 0| + -I - |4 >< 4| + |5 >< 5| + |6 >< 6| + |7 >< 7| + 



>< 



19 >< 91 



110 >< 101 



111 >< 111 



112 >< 121 



In (6), E is the unit matrix, and pa is the deviation matrix. The four states, \k >, 
k = 0,1,2,3 are the "active states", which are supposed to be manipulated without 
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noticeable change in the other twelve states. In binary notation, the "active states" 
are \00ij >. They are expected to evolve like the pure quantum states of a two-spin 
molecule. In Fig. 1, the "active states" are marked by a "bullet". 

The simplest CN gate for the system under consideration can be reahzed by a tt- 
pulse with frequency loq+SJ (see Fig. 1). For our case, we shall "invert" the convention. 
Namely, we shall associate with the ground state, |0 >, of the spin, the value "l"of the 
qubit; and the excited state, |1 >, with the value "0". 

A TT-pulse with the frequency {ujq + 3J) is expected to drive the right spin only if 
the neighboring spin is in the ground state, i.e., 

' >, iit = 0, (j=l-j), ^ ^ 

ignoring the overall phase factor. The frequency (u;o + 3 J) is unique frequency in this 
system. So, no other states are supposed to be changed by this pulse. 

To analyze the dynamics of the CN gate, we first calculate the 16 diagonal and 64 
non-zero off-diagonal matrix elements of the Hamiltonian (2). Next, we solve numer- 
ically the system of 256 equations of motion for the matrix elements of the density 
matrix Pa(^), for the initial conditions, 

^ 1 / 

E rnk{0)\n >< A;| + - - |4 >< 4| + |5 >< 5| + |6 >< 6|+ (8) 

n,k=0 ^ V 

|7 >< 7| + |8 >< 8| - |9 >< 9| - |10 >< 10| - |11 >< - |12 >< 12| . 

The coefficients r„,fc(0) (0 < n. A; < 3) in (8) describe arbitrary initial conditions for 
the "active part" of the deviation density matrix, pA- Certainly, 

3 

E^nn = l. (9) 
n=0 

The deviation density matrix, 

PA(t) = ^^^ E r^k{t)\n><k\, (10) 

^ n,k=0 

with the initial condition (8), was calculated numerically without any further approx- 
imations. 

The results of these numerical calculations are shown in Fig. 2 and 3, for the 
following values of parameters, 

ujk^l00{k + l), J =10, Q = 0.1, a; = 130, (11) 



Pa(0) 



^ Ea:=0 

2kBT 



4 



where /c = 0, 1, 2, 3. The characteristic dimensional parameters can be obtained, for 
example, by multiplying the parameters in (11) by 27rMHz. Fig. 2a corresponds to the 
initial conditions, 

roo(0) = l, r„fe(0) = 0, n,k<3, {n,k)^ (0,0). (12) 

One can see from Fig. 2a, that under the action of a 7r-pulse, the coefficient roo(t) 
(curve 1) decreases monotonically from 1 to 0, indicating the transition, 

|0000 >^ |0001 > . (13) 

Curves 2 and 3 in Fig. 2a show the evolution of the imaginary and the real parts of 
the average spin, Im < I'^ > and Re < I'^ > , correspondingly, where 

^+ = E4^ It -11 + ill. (14) 

0=0 

</+>=Tr{/+pA(t)}- E Itkrknit). 

n,k=0 

The value of Im < I'^ > describes the precessing amphtude of the spin which is 
maximum when the average spin is in the xy plane. 
Fig. 2b corresponds to the initial conditions, 

rii(0) = l, r„fe(0) = 0, n,A;<3, (n. A;) ^ (1, 1), (15) 

and describes the transition, 

|0001 >^ |0000 > . (16) 
Fig. 2c corresponds to the non-resonant initial conditions, 

r22(0) = l, r„fe(0) = 0, n,k<3, (n,k) ^ (2,2). (17) 

Fig. 2d corresponds to the initial conditions, 

r33(0) = l, r„fe(0) = 0, n,k<3, (n. A;) ^ (3, 3). (18) 

For both cases, shown in Fig. 2c and Fig. 2d, the population of the "active states" , 
\00ij >, does not change with accuracy of 10~^. The same is true for the population 
of all other "passive states" . 

Now we assume that under the action of a sequence of electromagnetic pulses we 
have a superpositional quantum state of a two-spin molecule, 

*(0) = coo(0)|00 > +coi(0)|01 > +cio(0)|10 > +cn(0)|ll >= (19) 
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co(0)|0 > +ci(0)|l > +C2(0)|2 > +C3(0)|3 > . 

For an ensemble of four-spin molecules at room temperature the superpositional state 
(19) corresponds to the initial values of r„fe(0) in (8), 

r„fc(0) = <(0)cfc(0), n,k<3. (20) 

Wc investigated the possibility of implementating a "tt"- pulse quantum CN gate for 
the ensemble of four-spin molecules. 

In Fig. 3 we show the evolution of the diagonal elements, rkk{t), under the action 
of a TT-pulse, for the initial conditions, 

co(0) = 0.3^/2^ c^{0) = 0.2'/^ C2(0) = 3-^/^ 03(0) = 6"^/^ (21) 

One can see from Fig. 3, that at the end of a 7r-pulse we have: roo = '"ii(O), rn = roo(O), 
and the values of and rss do not show a noticeable change. The same is true for 
the "non-active" diagonal elements, r44,..., ris 15. Thus, a "tt"- pulse provides the 
implementation of the quantum CN gate for the superpositional initial conditions (20), 
(21) as well as for the digital initial conditions (12), (15), (17), (18). 

We conclude that the dynamics of quantum CN gate can be separated into two parts 
- resonant and non-resonant, and the non-resonant interaction plays a minor role when 
the resonant condition u = ujq is satisfied, and the constant of interaction, J, is small 
enough. In this case, our numerical calculations confirm the possibility of implementa- 
tion of the quantum CN gate in an ensemble of quantum four-spin molecules, at room 
temperature, using a 7r-pulse. Further numerical investigations are required to analyze 
the dynamical behavior of the system as a function of the frequencies, Uk, in (11). 
When the frequencies Uk in (11) are close to each other, the resonant and non-resonant 
dynamics are expected to become more mixed. Also, additional numerical experiments 
are required to analyze the creation of the initial state (6), and to study the influence 
of sequence of 7r-pulses on the ensemble of spin molecules. These calculations are now 
in progress. 
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Figure Captions 



Fig. 1. The energy levels for a four-spin molecule. "Bullets" indicate the "active 
states" of the effective "pure" quantum two-spin system. 

Fig. 2. Evolution of the diagonal density matrix elements, r„„(t), in (10), and the 
average value, < >, for the parameters given in (11). (a) corresponds to the initial 
conditions (12); curve (1) refers to roo(t); (b) corresponds to the initial conditions (15); 
curve (1) refers to ru{t); (c) corresponds to the initial conditions (17); curve (1) refers 
to r22(t); (d) corresponds to the initial conditions (18); curve (1) refers to r^^lt). In 
Figs 2a-2d curve (2) refers to Im{< >), and the curve (3) refers to Re{< >). 
Vertical arrows indicate the end of the 7r-pulse. 

Fig. 3. Evolution of the diagonal density matrix elements, r„„(t), in (10) for the pa- 
rameters given in (11), and for the superpositional initial conditions (20)-(21). Curve 
(1) refers to roo(t); curve (2) refers to rii(t); curve (3) refers to r22(t); curve (4) refers 
to r33(t). Vertical arrow indicats the end of the 7r-pulse. 
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|1tH^"|14> 

|1(M1>^11> 



r11(HF-]1^ 



-|4> 



